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Classical and quantum phase transitions involve observables which are non-analytic as functions of a controlled 
thermodynamical variable. As occurs with the self-consistent Fermi Golden Rule, one condition to obtain the 
discontinuous behavior is the proper evaluation of a classical or quantum thermodynamic limit. We show that in 
presence of an environment, the oscillatory dynamics of a quantum two-level system, in analogy with a classical 
damped oscillator, can undergo a quantum dynamical phase transition to a non-oscillatory phase. This is obtained 
from a self-consistent solution of the Generalized Landauer Biittiker Equations, a simplified integral form of 
the Keldysh formalism. We argue that working at each side of the transition implies standing under different 
paradigms in the Kuhn's sense of the word. In consequence, paradigms incommensurability obtains a sound 
mathematical justification as a consequence of the non-analyticity of the observables. A strong case is made upon 
the need to deepen the public's intuition and understanding on the abrupt transition from static to dynamical 
friction regimes. 
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1. Introduction 

From metal and glass melting to steam en- 
gines, phase transitions have nurtured both tech- 
nological and scientific progress. Only in the 
last century has it become clear that phase tran- 
sitions occur when the relevant free energy is 
non-analytic on some controlled thermodynam- 
ical variables such as the temperature. How- 
ever, it is not trivial to see how this collective 
behavior emerges from the fundamental interac- 
tions governing microscopic variables. More re- 
cently, quantum phase transitions pQ, which are 
much more elusive, have received an increasing 
interest. This is mainly motivated by the High- 
Tc superconductors, transport in heavy fermion 
compounds and organic conductors and conden- 
sation of bosonic fluids. The greatest achievement 
was perhaps the theory of the extended-localized 
transition of electrons in disordered solids discov- 
ered by P. W. Anderson in 1958 [2]. The whole 
solid state community was taken by surprise by 
his statement that non-interacting electronic or 
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vibrational eigenstates in solids would transform 
from Bloch plane waves into exponentially local- 
ized functions whenever the strength of a homo- 
geneous disorder exceeds a critical value. Typi- 
cally, one should recognize a phase transition as a 
non-analytic behavior of the ground state energy 
or other observable as a function of a control pa- 
rameter g [jQ . This sort of phenomena could triv- 
ially occur when the control parameter moves the 
system through a level crossing. However, this in- 
volves a total Hamiltonian of the form Hi + gH 2 
with Hi and H 2 mutually commuting. In finite 
systems this would be an extremely rare situa- 
tion, but it becomes more likely when one con- 
siders an infinite lattice. In this case, the infi- 
nite number of degrees of freedom involved could 
transform an avoided crossing of the finite system 
into an actual level-crossing. 

In this work we want to discuss how a quantum 
dynamics of a system can undergo a phase transi- 
tion. We consider a system tunneling coherently 
between two levels to form a Rabi oscillation. 
This system is ubiquitous in Nature [3J, but has 
received renewed attention in quantum informa- 



2 



tion field because it constitutes a swapping gate 
4 5 6J. The presence of a quantum environment, 
requires the solution of the dynamics of open sys- 
tems [7|- We resort to the Keldysh formalism [5] 
which, with some simplifying assumptions, be- 
comes the Generalized Landauer-Biittiker Equa- 
tions |9|10j which can be solved analytically. We 
find that the oscillatory dynamics can freeze when 
the interaction with a quantum environment ex- 
ceeds certain critical strength. This behavior has 
a close analogy with the transition between dy- 
namical regimes (oscillating-overdamped) under- 
gone by a classical oscillator when friction is in- 
creased. Since several of the current descriptions 
of these phenomena do not point out the concep- 
tual assumptions enabling the phase transition, in 
this article I will sketch out the calculations focus- 
ing on the conceptual conundrums: What is the 
meaning of a 'thermodynamic limit' in classical 
and quantum mechanics? Why does the quan- 
tum description of an open system involve a form 
of thermodynamic limit, and why can this enable 
a quantum dynamical phase transition? 

Finally, I will conclude with a section associat- 
ing phase transitions to a paradigm shift in sci- 
ence [111 . Similarly to what occurred with the 
Aristotelian-Newtonian shift, the mechanicists- 
probabilistic shift manifested in the well known 
Loschmidt vs. Boltzmann polemics (that 
switches between reversible and irreversible me- 
chanics) and the related Zermelo/Poincare vs. 
Boltzmann argument on the transition between 
recurrent and dissipative mechanics |12j. 



2. Effective Hamiltonians 

We are particularly interested in the coherent 
polarization transfer among two magnetic nuclei, 
which can be reduced to a non-interacting elec- 
tron |13j so we will resume the basic formulation 
of the latter problem [141 . The real symmetric 
Hamiltonian H — + V, describes the dynam- 
ics of two states, \A) = c\\$) and \B) = |0) 
which are mixed by a tunneling matrix clement 
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Eliminating one of the amplitudes, e.g. ub, gives 
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Obviously the bracket is an "effective" Hamilto- 
nian H A S - — Ea(s) which includes the "energy 
shift" due to the eliminated orbital 



E A (e) = E A + E A (£), 



1 



-V, 



BA- 
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(4) 



^ e) = V **{e-E B ) 

Indeed, under an apparent simplicity, the equa- 
tion becomes non-linear and it's solution provides 
the two exact eigenvalues of the system 



e A = \[{Ea + E B ) - hw AB ], 

SB = \[{E A + E B ) + tuJAB] 

flUJAB = J{E B - E A ) 2 +4\Vab\ 



(5) 
(6) 

(7) 



This procedure can also be expressed in terms of 
Green's functions. Given a positive ry, one defines 
the retarded and advanced resolvent matrices, 



G R (s + u 1 )=[G A (e-ir,)] 1 

= [(e + ir))I -Hf 1 



(8) 
(9) 



(i - E A ) (e + ir)- E B ) - VabVba 



(e+irt-EA){e+iv-£B) 

e + ii]-E B -Vab 



-Vba 



£-E a 



(10) 



The retarded (advanced) Green's functions are 
matrix elements which, for real e, have diver- 
gences at the eigen- energies as rj — > + being ana- 
lytic in the upper (lower) half plane. These diver- 
gencies weigh the probability of the unperturbed 
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state on the eigenstates \A} and \B} . Hence, the 
local density of states (LDoS) at site n = A,B 
results: 



NJe) = -i lim Im (n\ G oR (e + m) In) (11) 

= \(n\ A)\ 2 S(e -e A ) + \ (n\ B)\ 2 S(e - e B ). 
The diagonal matrix elements can be rewritten 



as 



E A (e) 



(12) 



Identifying the unperturbed Green's functions 
G°^(e) = [e — En}^ 1 and expanding one gets, 



Gl A (s) 



1 



Gf A {e) 



X A (e) 



= Gf A (e) + G a r a {e)H a {e)G^ a {e) 
+ G° a r a (e)Z a (e)G a r a (e) 

xX a (e)G° a r a (e) + .... 



(13) 



This shows that the exact solution is the sum of 
an infinite geometric series. This is represented 
as Feynman diagrams in Fig. 1. There is yet 
another form of writing this, a Dyson equation, 



GXAs) = G AA (e) + G% a (e)V a {e)G a h a {e) 



(14) 



It is useful to note that all the above results, 
and most of what follows, are also valid if \A) 
and \B) denote whole subspaces. In that case, 
all the presented equations and diagrams hold 
but with matrix elements transformed themselves 
into matrices |15|16j . We might choose not to deal 
explicitly with an undesired subspace, for exam- 
ple the whole subspace \B) , and still get and ef- 
fective Hamiltonian restricted to the subspace \A) 
and also the exact Green's function. 

Usually, given an initial state, the dynamics 
is evaluated from eigen-energies and eigenstates. 
Alternatively, it can be expressed in terms of 
Green's functions. For example, the probability 
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Figure 1. In the upper panel, the exact Green's 
function ( thick line) is represented as an infi- 
nite series of unperturbed Green's functions ( thin 
lines). Coupling matrix elements are dashed lines. 
The botton panel shows the self-consistent Dyson 
equation and the self-energy. 



that a particle which was in the state \A) at t = 
is found at state \B) at a later time t results: 



P BA {t) = (Blexpl-Lff t]\A) 0[t] 



(15) 



7^0+ J Ann 



= J d£P B ^ A (e,t), 



(16) 
(17) 

(18) 



with 



P B , A (u>) = J dE^G R >a (e + \?vu>)G% b {e - ±M- 

(19) 

The appearance of the function 9[t] in Eq. [15] 
is consequence of the election of the sign of the 
imaginary part in the retarded Green's function. 
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Figure 2. (Color online) Energy-time distribu- 
tion function for a two-level system (in units of 
V and h/V respecively) The dark (yellow-red on- 
line) and clear (blue online) regions differ in sign. 
The formed stripes manifest the progressive de- 
crease in the small structure's scale as function of 
time. 



The remaining two lines constitute alternatives 
for writing the product of the independent in- 
tegrals. The function PB,A{s,i) (as well as its 
transform Pb,a(s,u)) is not an actual probabil- 
ity but a form of energy-time distribution func- 
tion from which a real probability can be obtained 
as a marginal distribution, i.e. by integration 
of one of the variables. In more general prob- 
lems, this energy-time distribution enabled [9110] 
to consider time dependent statistical distribu- 
tion functions. For the particular case of equal 
energies Ea = Eb = and Vab — V with the su- 
perposition | A) = ^7= (| A) + \B)) as initial state: 



fU,A(e,t) 



sin [2 (e + V) t] 6(t) (20) 



V 2 +e(V- 2s) 
2e{V^-t 

V le(v^) sm[2{e-V)t]6{t) 



This distribution oscillates as a function of each 
participant energy at a rate which is determined 
by its distance to the eigenvalue (see Fig. 2). 
From it, the Rabi oscillation is obtained as: 



Notice that while the result of the integral re- 
mains a simple oscillation, for long times the in- 
tegrand becomes an increasingly oscillatory func- 
tion on the energy variable. In a numerical inte- 
gration, regions with too small structures would 
contribute as pseudo-random amount to the in- 
tegral making it numerically unstable. It would 
be tempting to do an analogy with similar struc- 
tures in the standard momentum-position Wigner 
function suggested by Zurek [17118] . and inter- 
pret this phenomenon as a manifestation of the 
instability of this quantum superposition towards 
decoherence. In fact, ideal Rabi oscillations con- 
trast with experimental observations, such as Fig. 
4-a of Ref. [TH] , where the environment is actually 
attenuating the oscillation while the probability 
is conserved. Thus, our simple quantum mechan- 
ical model should be extended to include some 
form of environmental interaction. 



3. The spectrum of a finite linear chain and 
continued fractions. 

We will represent the environment with our fa- 
vorite model, the linear chain. It not only rep- 
resents a chain of spins interacting through a 
XY interaction [20] but it is a reasonable model 
for polymers, quasi-one dimensional crystals and 
metal wires. Even a crystal structure can be re- 
duced to a set of uncoupled linear chains. We 
start by adding a third state to our two state sys- 
tem, 



H 



Ei 


-V12 





-V 2 i 


E 2 


-v 23 





-v 32 


E?, 



(22) 



We start with V12 — 0. Through the identifica- 
tion of the indices 2 — * A, and 3 — * B, we use 
Eq]3] eliminate state B{ i.e. 3) so that G A A (e) 
— ^ Gvj^e). Now we turn-on Vi,2 and identify 



Pa A*) 



dsP A>A (s, t) = cos 2 auj A Bt). (21) 



1 -> A and 2 
of B to get: 

GUe) 



B, and we repeat the elimination 



(23) 



[Gl R M -Vi 2 G° 2 *(e)V, 
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We replace it and obtain a nested fraction: 



4. The semi-infinite ordered chain 



e-Ei- V r . 



-V21 



e — E 2 — V23 - 



:\2 



S 2 



(24) 

In the present context, the self-energy accounts 
for presence of states at the right. 

Hamiltonian of Eq. ([22]) presents an interesting 
phenomenon. If V23 -C V12 = Vab the system AB 
is well defined and site 3 can be seen as an "envi- 
ronment" weakly perturbing the system through 
Vse = V23. If we allow the parameters to switch 
to the opposite regime Vse = V23 > V12 = Vab, 
state B becomes "captured" by the environment 
and the state A becomes almost isolated. This 
can be seen as a form of the Quantum Zeno Ef- 
fect [2T] caused by the internal degrees of freedom 
[22123] , 

Since the procedure performed above was 
in fact a step of a renormalization group 

the 



algorithm [24(15], we can iterate it to get 
general continued-fraction that describes a chain 
with N orbitals: 



En ( £ ) ^ Vn,n+1 ' 



1 



E n — £ n +i (e) 
together with the termination condition 

Sat (e) = 0. 



Vn+\,n- (25) 



(26) 



Hence, the Green's function, as the self-energy, 
is the ratio between two polynomials. This yields 
the N eigen-energies and eigenvalue weights of 
the finite system. As predicted by Poincare this 
produces many recurrences. A particularly inter- 
esting dynamical recurrence is what we called [25] 
the mesoscopic echo which appears at the Heisen- 
berg's time Tme = h/ A where A is the mean 
level spacing. Signatures of this phenomenon 
where experimentally observed in Cordoba 20J 
and confirmed in Zurich by the group of Richard 
R. Ernst as can be seen in Fig. 3-B of Ref . [55] . 



When the chain of lattice spacing a is ordered 
(E n = 0, V n>n -\-i = V) and infinite there is no 
termination condition as Eq. [221 Instead, all sites 
"see" the same environment at their right. Hence, 
the equation that is now equivalent to the Bloch 
theorem is 



£„ (e) = £„+i (e) = £ (e) , 
from which: 



E(e) 



V 2 



(27) 



(28) 



£-£(£) 

The surprise is that in the region where there are 
real eigenvalues, the solution is complex 

E( e )=A(e)-iT( e ), (29) 

the energy shift is a piece-like function: 



A(e) 



(I)' 



V 2 for e > 2 \V\ , 
for |e|<+2|V r |, 



-J - V 2 for e < -2 |V| . 



(30) 



while the group velocity, T — hv e /a, results 







r( e ) 



V 2 




(1)' 



for e > 2 \V\ , 

for |£|<+2|y| 
for £ < -2 |^| . 



(31) 



The sign of the square root is consistent with 
the analytical properties described above, while 
the real part goes to zero as lim^ioo A(e) = 
which means that the spectrum of the linear chain 
remains bounded after the interaction has been 
turncd-on. The consistency of these solutions can 
be checked through the convergence of the self- 
energies in chains of increasing lengths. This ex- 
presses the Quantum Thermodynamic Limit: 



-T (e) = lira lim ImE 1 (£ + i?7) 
^ lim lim ImS, (e + iri) 



(32) 
(33) 
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a.e.e means for almost every s, i.e. except for a 
set whose probability measure is zero. The non- 
triviality of this limit is manifested in the fact 
that it is non-uniform. 

5. The Fermi Golden Rule as a Quantum 
Thermodynamic Limit 

In the above discussion we obtained an effective 
energy with an imaginary component. It actu- 
ally means that perturbation theory does not con- 
verge. The unperturbed eigenstate is so far from 
the new eigenstates that their scalar product van- 
ishes. In the dynamics, this should manifest as a 
progressive decay where the Poincare recurrences 
no longer appear. This means that the prob- 
ability escapes towards the semi-infinite chain. 
For the homogeneous linear chain this involves 
a power law decay according to the law Pi,i(t) — 
(Vt) . A particularly interesting case occurs 
when at the end (surface) of this semi-infinite 
chain we add an orbital (or atom) with energy 
E and interaction Vq <C V. This adatom model, 
is a particular case of the Friedrichs model. One 
knows that this situation leads to a typical expo- 
nential decay described by the Fermi Golden Rule 
(FGR). However, a deeper analysis shows that the 
exact rate of decay differs from that in the FGR. 
The new rate, Tq/H, arises from a Self Consistent 
Fermi Golden Rule [27] ■ It is the imaginary part 
at the exact pole e r — iT of the Green's function: 



than, 



ir 



V 2 



ir„ 



(34) 



which can be obtained analytically or by itera- 
tion. 

One should not forget that a quantum decay 
starts always quadratically, in this case with a 
time scale H/Vq. It only starts looking exponential 
after a time t$- This is a short time scale, 



t s = hir Ni(e r ), 



(35) 



when the escape from the surface site towards the 
rest of the chain prevents the return and hence 
stops giving an appreciable contribution to the 
survival. Here, Ni(e r ) is the LDoS at the surface 
site in absence of the adatom. At times longer 



t R 



a— In 

n 



4 

1 



(36) 



the return amplitude, determined by the high or- 
der processes that has already escaped but re- 
mains in the neighborhood, starts being compara- 
ble to the pure survival amplitude. From then on, 



-3/2 



Here. 



decay becomes a power law [r(e r )i] 
B = 4V is the bandwidth and a, (3 > 1 are con- 
stants that depend on the van Hove singularities 
of N±(e r ) and other details of the model. At tn a 
striking destructive interference between the pure 
survival amplitude and the return amplitude may 
occur. In quantum systems, this "survival col- 
lapse" [27] has yet to be observed. 

In summary, the validity of the FGR is re- 
stricted by memory effects to times between tn 
and ts- The standard FGR holds in the wide 
band limit Tq/B — > which also implies that 
VoNi(e r ) — > 0. It is only in this condition, valid 
in a quite broad variety of situations, that one 
can forget the quantum memory effects of a fi- 
nite bandwidth and replace both A (e) — ir (e) 
by A — ir independent of s. The environment 
behaves as a Markovian process and we refer 
to them as the "broad band approximation" or 
"fast fluctuations approximation" . One should 
be careful, however, interpreting this as an "ir- 
reversible" behavior [28] . Actual irreversibility is 
consequence of an instability that manifests when 
one attempts to revert the dynamics through a 
"Loschmidt daemon" , H — > —H [29]. One should 
note that an imaginary part by itself does not 
ensures irreversibility as long a one can change 
the sign of the whole Hamiltonian. The instabil- 
ity of this procedure can be tested and quanti- 
fied through the Loschmidt echo (or fidelity) in 
spin dynamics [30 31 , confined atoms [32] and 
microwaves in cavities [33] • See also Ref. [34] for 
a completely different approach to achieve time- 
reversal. 

The physical meaning of the imaginary part we 
introduced at the beginning is now evident: it 
represents the weak interaction with an environ- 
ment. In such situation, A does not contribute 
much to the dependence on e and one includes 
it by shifting the energies. This approximation 
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would give a steady decay of the Rabi oscillation 
as is indeed typical of many experimental setups. 
See for example the Fig. 4-a in Ref [TIJ] • However, 
one might wonder how to return the probability 
lost in this decay. In fact in presence of two iden- 
tical linear chains connected to states A and B, 
one would find probability [35j . 

P A>A (t) = P A>A (t) exp [-2Tt/K] (37) 
= cos 2 (2ujABt) exp [— t/r] , with r = h/2T. 

Clearly, this describes the evolution of polariza- 
tion tunneling between two nuclei shown in Fig. 
4.b of Ref. 19 1 . In this case, the probability 
(polarization) is not conserved but it decays ac- 
cording to the FGR. While this could be correct 
in some physical situations, the description of a 
situation closer to Fig. 4. a, where probability is 
conserved, remained a challenge. 

6. The Generalized Landauer-Biittiker 
Equation 

The imaginary energy has been a puzzle for ev- 
eryone using Green's functions and regularizing 
its poles. Sometimes, as in the electron-phonon 
processes, an explicit form for this imaginary en- 
ergy is evaluated through the FGR. Even the 
transport equations, as the Kubo formula, rely 
on some natural broadening which enables the 
computation but produces local non-conservation 
of currents. The answer was given by D'Amato 
and Pastawski [36] who, extending an idea of 
Buttiker[37], realized that the escape to an envi- 
ronment is equivalent to saying that, at each time, 
a fraction of the system occupation escapes to the 
chain which could act as a voltmeter. As an ac- 
tual voltmeter, however it should not extract net 
particles from the system, so it returns a particle 
for each one collected. This can be expressed [9] 
in terms of the Landauer description of transport 
which now accounts for time dependences and 
decoherent process in the form of a Generalized 
Landauer-Biittiker Equation (GLBE). Hence, for 
every process of "escape" from the coherent beam 
due to the interaction with the environment, a 
fresh incoherent particle must be reinjected into 
the system as expressed in Eq. (3.7) of Ref. [9]. 



-<<X> = ^<d"> + <^jfc> 

n 

m n 

+ ^<U> B + ••• 

m n 

n 
n 



Figure 3. Diagrams for the density propaga- 
tor from A to B as dictated by Generalized 
Landauer-Biittiker Equation. Horizontal lines are 
single particle Green's functions dressed by the 
environment. Shadowed vertices are the self- 
consistent density propagators. The vertical dou- 
ble dashed lines represent the reinjection pro- 
cesses. The last collision occurs at site n. 



This physical picture, finds its formal justification 
when the system-environment interactions are lo- 
cal and the environment spectrum is so broad 
that it becomes instantaneous and energy inde- 
pendent. In this case, the Keldysh quantum field 
theory formalism, expressed in its integral form 
[38] , reduces to the GLBE [10114] represented in 
Fig. 3. 

We consider a degenerate two-level system 
where, besides tunneling, each state suffers the 
action of a complex self-energy, such as that of 
the linear chain described in the previous section. 
This results in an homogeneous interaction with 
the environment. The natural parameter regulat- 
ing the effectiveness of the system-environment is 
g = uj A bV jh. For this two-level system the GLBE 
becomes: 

?A,A(t) = p A,A( t ) (38) 
rU _ dt ~ 

+ V / P A ,n{t-ti)— P n , A {ti), 

and equivalent equations for the components BA, 
AB and BB. This is a Volterra's type equation. 
This is a Dyson equation (much as Eq. [T4")) for 
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a density, i.e. a two-particle Green's function, 
and is also known as a Bethe-Salpeter equation. 
The first term describes the probability of coher- 
ent propagation from the initial to the final state 
which decays due to interactions with the envi- 
ronment. The kernel of this equation is precisely 
PA,n(t — U), the two-particle propagator. Since 
dti/r is the probability of having the last inter- 
action with the environment at the time inter- 
val dti around ti. The solution of the homoge- 
neous GLBE can be obtained by Fourier trans- 
formation [39] Pb,a{w) and decays toward the 
equilibrium Pb,a{^) — > \- One notable thing 
is that the first term in the right has poles in 
the complex w-plane that correspond to the dif- 
ference of energies and do not present any form 
of non- analyticity. The self consistent solution 
Pb,a{u) has more information. In fact, the poles 
of SPb,a(u) = Pb,a{w) — \5{uj) are precisely at 

^ ± -ir = ± v / [w AB ] 2 -r 2 -ir (39) 

The trajectories in the complex plane are shown 
in Fig. 4-b. The important feature is that the 
real part of the poles (Fig. 5-a) collapses at for 
a critical value g c — 1 and from this point they 
split in two terms of null real part. One of them 
decreases with environment interaction whereas 
the other decreases. It is the later that controls 
the long time behavior. 

SP AA (t) = P A At) - 3 = a o cos [(w + ir )* + 4>] ■ 

(40) 

Here Pj±\ = \ is the equilibrium occupation 
while a\ = [4uj 2 t 2 + l] / (16cj 2 t 2 ) and </» = 
arctan [l/2wr] warrant the initial cuadratic de- 
cay. 

The described behavior has been experimen- 
tally observed in heterogeneous polarization 
transfer, see Fig. 7 of Ref. [30], but overlooked 
since the early theory for this experiment [40] did 
not contain the transition. Recent experimen- 
tal studies [5] show the divergence of the period 
2it/oj at a critical ratio Huiab/^ = 1- Concep- 
tually, the transition is from an isolated system 
that is weakly perturbed to a state in which the 



a) G*(e) b) 5P AA (a) 




Figure 4. a) Paths of poles of a single particle 
Green's function, e.g. G AA (s), when an homoge- 
neous decay Tis increased. They move parallel to 
the imaginary axis, b) Paths of poles of the ob- 
servable 5Paa(&) (a two-particle self-consistent 
Green's funcion) when T increases. The sym- 
metric frequencies collapse at the center where a 
branching occurs. One mode becames long life 
while the other has a short life time. 



effect of the environment is no longer perturba- 
tive. The system may be no longer well defined as 
discussed with reference to Eq. [22j This would be 
a dynamical Quantum Zeno Effect [22123] , While 
the limiting cases were somehow expected, it was 
by no means obvious that this change could be 
critical. The non-analyticity was enabled by the 
infinite degrees of freedom of the environment 
in the proper quantum thermodynamic limit and 
the self-consistent nature of Eq. [33 

The excess density 8PA,A{t) behaves exactly as 
the amplitude x(t) of a damped classical oscil- 
lator which undergoes a transition to an over- 
damped regime. Indeed, considering a damped 
harmonic oscillator of mass m and natural fre- 
quency uj , the relaxation rate as a function of the 
friction coefficient, T, follows precisely the trace 
of Fig 5-b: The rate increases with friction un- 
til a critical value when it starts to decrease in- 
versely proportional to the friction coefficient T. 
This is, of course, a non-analytic critical behav- 
ior. Where does it come from? From the imag- 
inary self-energy correction that shifts the nat- 



ural frequency uj in the oscillator's dynamical 
susceptibility x( w ) = —m^ 1 / [uj 2 — (iv 2 — ic^r)] . 
The damped Newton's equation is not a funda- 
mental law but it is written on phenomenolog- 
ical grounds. However, the inclusion of T can 
be justified, within statistical mechanics, by in- 
cluding the action of a Brownian bath [41]. Re- 
cently, we obtained a simpler demonstration 42J 
using as environment a chain of oscillators whose 
N degrees of freedom are considered by taking 
the thermodynamic limit of N — > oo precisely 
in the same way as described above in the con- 
text of the FGR. It is interesting to note that 
while 2cj Q /r ^> 1 corresponds to the standard os- 
cillation. In a similar way, in the quantum case 
2w Q /r 3> 1, the system is well defined and the 
environment is a small perturbation. In contrast 
in the regime controlled by friction, 2lu /T <C 1, 
the inertia term can be completely neglected. 

It is clear that most of the qualitative fea- 
tures of the spectral properties described above 
are valid for other linear systems (provided that 
there is a thermodynamic limit) and hence are 
ubiquitous in Nature. In magnetic resonance, a 
phenomenon known as exchange narrowing, has 
long been described and clearly observed [H] . 
However, its explanation requires either Brown- 
ian fluctuations or the use of Markov chains with 
imaginary probabilities. . . ! [45] . 

7. Phase Transitions as Paradigm Shifts 

In the previous sections we have touched upon 
issues such as complex energies, imaginary prob- 
abilities, irreversibility, recurrences, decoherence, 
non-analytic observables, etc., all of them gen- 
erating strong polemics. In consequence, some 
epistemological comments are pertinent. 

One of the central statements of ancient 
Physics was Aristotle's dictum that everything 
that moves is moved by something else. More pre- 
cisely, Aristotle says that the velocity of a mov- 
ing object is directly proportional to the force 
and inversely proportional to the resistance, i.e. 
x = F/T. In the absence of a proximate force, 
the body would come to rest immediately. Obvi- 
ously, a difficulty found in the Aristotelian view is 
the justification of why a projectile keeps moving 
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Figure 5. a) The frequency of a two-level system 
(Fig. 4b) collapses at zero for a critical T. b) the 
decoherence rate as function of T. This also rep- 
resents a relaxation rate in a damped harmonic 
oscillator as function of friction strength. The 
botton left point is ideal frictionless Hamiltonian 
mechanics or Newton's paradigm. The right side 
is the realm of Aritotle's paradigm where inertia 
becomes negligible. 



through the air. The logic of the explanation is 
not as clean as the central statement: a projectile 
would owe its continuing motion to the force of 
eddies or vibrations in the surrounding medium, 
a phenomenon known as antiperistasis. This was 
formalized later on by the scholastics [47] who 
proposed that motion was maintained by some 
property of the body, the impetus, which once 
set in motion, would impart the force keeping the 
movement. Buridan's impetus has the same con- 
sequence, but very different justification, than the 
modern concept of momentum |46] . 

Physics seemed a quite solid construction un- 
til the experiments and intuition of Galileo and 
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analytical skills of Newton started to build much 
of our current knowledge. In this new formula- 
tion, the inertia is no longer a correction but the 
fundamental principle. Friction appears now as 
a subsidiary and phenomenological effect needed 
to fit the natural phenomena to the ideal scheme. 
Indeed its justification remained for a long time 
alien to Hamiltonian mechanics. One had to 
wait for the appearance of Boltzmann's statisti- 
cal mechanics and the work of Smoluchowski and 
Einstein to have a place in the theory building 
(for a simple Hamiltonian model justifying fric- 
tion see Ref. [42] ). In any case, Aristotelian 
and Newtonian views, were so completely irrec- 
oncilable that Thomas Kuhn concluded that 
they were indeed different views of Nature. He 
coined the term paradigm shift to describe a 
change in the basic assumptions within the rul- 
ing theory of science. According to Kuhn, sci- 
ence progress relies not only on a linear accu- 
mulation of new knowledge as sustained by Karl 
Popper but, more fundamentally, on periodic rev- 
olutions in which the nature of scientific inquiry 
within a particular field is abruptly transformed 
[llj . Rival paradigms are said to be incommen- 
surable because it is not possible to understand 
one paradigm through the conceptual framework 
and terminology of another rival paradigm. 

What seems disturbing to some scientists, is 
the possibility that no ultimate truth underlies 
this confrontation between paradigms 49J . Is 
it possible to synthesize these extreme behav- 
iors into a single framework? Our answer is yes, 
because incommensurability involves comparing 
states at different sides of a phase transition. In- 
deed, Aristotle's paradigm is placed at the ex- 
treme right side of Fig. 5-b where the inertia's 
contribution to an equation of motion is com- 
pletely neglected. The impetus corrections allows 
one to move somehow to the left. The contrast- 
ing Newton's paradigm, x — F/m, is placed at 
the extreme left, in the origin of Fig. 5-b. From 
that ideal point one could conceive adding fric- 
tion as a correction. Consider a mass placed in a 
bowl where friction can be varied. Would anyone 
experimenting in one of those extremes conceive, 
without completing the experiment of the other 
regime, what the response at the other edge would 



be? The answer is a clear no, as the non-analytic 
function does not allow a natural extrapolation. 
Indeed, it was not until Gauss popularized the 
concept and interpretation of Euler's complex 
numbers that both regimes fitted into a single de- 
scription. Even with that tool, numerous discus- 
sions with students and colleagues convinced me 
that intuition fails lamentably at the non-analytic 
point. The same occurs when one discusses prob- 
lems which involve the non-homogeneity of the 
limits, which indeed is at the root of the micro- 
scopic description of friction. Many other con- 
troversies in Physics have a resolution within this 
framework: we have already advocated that the 
Loschmidt vs. Boltzmann controversy is a con- 
sequence of the non-uniformity of the limits for 
an imperfect time reversal experiment |31j . Each 
argument results valid in a different approach to 
the limiting case (see Fig. 6 in Ref. [H]). The 
Zermelo/Poincare-Boltzmann controversy is an- 
other consequence of different forms of taking the 
thermodynamic limit. 

More recently, in the quantum framework, the 
localized-extended transition owes its origin to 
the fact that strong disorder induces a non- 
uniformity of the limits respect to ensemble aver- 
age, 

f ( E ) = lim ( lim Im E(e + in)) cns (41) 

avc. 

7^ ( lim lim ImSfe-fi^) ) =0. 

\ 77—^0+ N— >oo / cns. a.c.e 

ave. 

(42) 

This inequality and the last equality were 
proved and tested numerically in Ref. [50] . They 
show that in the localized regime the spectrum is 
pure-point. Not recognizing it led to contradic- 
tory results for about two decades |51|52j . Also 
the coarse grain average has subtle properties 
of non-uniformity respect to the thermodynamic 
limit which need further exploration [53]. 

There are other smaller paradigm shifts in con- 
densed matter physics, which resulted somehow 
less confiictive, produced by the need to explain 
quantum phase transitions. We can mention 
superconductivity (from current carried by sin- 
gle electrons to Cooper's pairs), localization and 
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mesoscopic transport (which shifted from Kubo's 
view where dissipation occurs inside the sample 
to that of Landauer, where it occurs at the ex- 
ternal reservoirs) and the Integer Quantum Hall 
(where the standard vision of bulk current yields 
to Biittiker's edge current). 

Finally, I feel the obligation to mention an- 
other phase transition which should not be much 
different from that discussed above: the transi- 
tion from static friction to dynamical friction. In 
that case, ordinates in Fig. 5-b describe the fric- 
tion force as a function of the applied force. The 
abrupt fall of the last at a critical force describes 
the transition to the almost constant value of the 
dynamical friction. In fact, the non-analytic jump 
from static friction to dynamical friction is so 
unexpected and counter-intuitive that no other 
phase transition seems to have a bigger deathly 
tall in "accidents" on the road, at work or even 
at home. It seems to me that it is a most urgent 
challenge to devise an educational strategy capa- 
ble to develop, in the general public and physi- 
cists alike, an intuition on this phenomenon. On 
the physical side, friction has only recently been 
reintroduced as a fundamental problem [54] . Its 
formulation relies on models having a close con- 
nection to issues discussed above [55] . This is still 
another phase transition that opens new ques- 
tions not only for basic physics but, even more 
importantly, also to social and cognitive sciences. 
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